ID multicomponent Fermions with delta function interaction in strong and weak 

coupling limits: K-component Fermi gas 
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We derive the first few terms of the asymptotic expansion of the Fredholm equations for one- 
dimensional K-component fermions with repulsive and with attractive delta-function interaction in 
strong and weak coupling regimes. We thus obtain a highly accurate result for the ground state 
energy of a multicomponent Fermi gas with polarization for these regimes. This result provides 
a unified description of the ground state properties of the Fermi gas with higher spin symmetries. 
However, in contrast to the two-component Fermi gas, there does not exist a mapping that can unify 
the two sets of Fredholm equations as the interacting strength vanishes. Moreover, we find that 
the local pair correlation functions shed light on the quantum statistic effects of the K-component 
interacting fermions. For the balanced spin case with repulsive interaction, the ground state energy 
obtained confirms Yang and You's result [Chin. Phys. Lett. 28, 020503 (2011)] that the energy per 
particle as k —¥ oo is the same as for spinless Bosons. 
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I. INTRODUCTION 

The recently established experimental control over the 
effective spin-spin interaction between atoms is open- 
ing up new avenues for studying spin effects in low- 
dimensional atomic quantum gases with higher spin sym- 
metries. Fermionic alkaline-earth atoms display an ex- 
act SU(k) spin symmetry with k = 21+1 where I 
is the nuclear spin For example, a recent experi- 
ment Q for m Yb dramatically realised the model of 
fermionic atoms with SU (2) <25 SU (6) symmetry where 
electron spin decouples from its nuclear spin / = 5/2. 
Such fermionic systems with enlarged SU(k) spin sym- 
metry are expected to display a remarkable diversity of 
new quantum phases and quantum critical phenomena 
due to the rich linear and nonlinear Zeeman effects. The 
study of low-dimensional cold atomic Fermi gases with 
higher pseudo-spin symmetries has become a new fron- 
tier in cold atom physics. 

One-dimensional (ID) quantum Fermi gases with 
delta-function interaction are important exactly solvable 
quantum many-body systems and have had tremendous 
impact in quantum statistical mechanics. The spin-1/2 
Fermi gas with arbitrary polarization was solved long 
ago by Yang Q using the Bethe ansatz (BA) hypoth- 
esis. Sutherland 3 generalized the result of the spin-1/2 
Fermi gas to ID multi-component Fermi gas in 1968. The 
study of multicomponent attractive Fermi gases was ini- 
tiated by Yang [3 and by Takahashi @ . Using Yang and 
Yang's method Q for the boson case, Takahashi Q and 
Lai |9J derived the thermodynamic Bethe ansatz (TBA) 
equations for spin-1/2 fermions. In the same fashion, 
Schlottmann [1J| derived the TBA equations for SU(n) 
fermions with repulsive and attractive interactions re- 
spectively. Recently, the thermodynamics of the multi- 
component Fermi gas was obtained by solving the TBA 



equations in [l2j. These models with enlarged spin sym- 
metries have received a renewed interest in cold atom 
physics [lltlMII. 

In this communication, we consider ID K-component 
fermions with repulsive and with attractive delta- 
function interactions. Although the model was solved 
long ago by Sutherland Q, solutions of the Fredholm 
equations for the model are far from being thoroughly in- 
vestigated except for the two-component Fermi gas [i~6| . 
From a theoretical point of view, finding a general form 
of the solutions to the Fredholm equations for multicom- 
ponent Fermi gases with higher spin symmetry imposes 
a number of challenges. 

In the present paper, using the method which we re- 
cently developed in the analytical study of the ID two- 
component Fermi gas [la ], we approximately solve the 
Fredholm equations of the ID /{-component fermions 
with polarization for the A) strongly repulsive regime; 
B) weakly repulsive regime; C) weakly attractive regime; 
and D) strongly attractive regime. We thus obtain the 
ground state energy of the K-component Fermi gas with 
polarization that provides a fundamental understanding 
of the ground state properties, such as phase diagram, 
magnetism and quantum statistical effects. In contrast 
to the two-component Fermi gas [l6j], the two sets of 
Fredholm equations for the multi-component Fermi gas 
with weakly repulsive and attractive interactions can not 
be unified by the density mapping as the interacting 
strength vanishes. The multiple spin degrees of free- 
dom impose subtle intricacies of quantum statistics in the 
ground state properties of the systems. We further study 
the local pair correlation functions for the K-component 
interacting fermions. We find that the local pair correla- 
tion as k — y oo is the same as for spinless bosons. This 
result is consistent with Yang and You's finding [l9| that 
the energy per particle as k — > oo is the same as for 
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spinless Bosons. 



the ground state, the generalized Fredholm equations for 
c > are given by 



II. THE FREDHOLM EQUATIONS 



The Hamiltonian for the ID iV-body problem is [H, H[ 



H = - 
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(1) 



It describes N fermions of the same mass m confined to a 
ID system of length L interacting via a 5-function poten- 
tial. There are n possible hyperfine states |1), |2), . . . ,\k) 
that the fermions can occupy. For an irreducible rep- 
resentation [k n ",(k- l)^- 1 ,.--^ 2 ,^ 1 ], the Young 
diagram has k columns with the quantum numbers JVj = 
N l — N 1+1 , here the iV 1 is the numbers of fermions at 
the i-th hyperfine levels such that N 1 > N 2 > . . . > N K . 
This system has SU(n) spin symmetry and U(l) charge 
symmetry. The coupling constant gm can be expressed 
in terms of the interaction strength c = —2/a±D as 
9id = h 2 c/m where a\u is the effective ID scattering 
length. Here c > for repulsive fermions, and c < for 
attractive fermions. 

The energy eigenspectrum is given in terms of the 
quasimomenta \k{\ of the fermions via E = X^jLi kj> 
which in terms of the function eb(x) = (x + ibc/2)/(x 
ifec/2) satisfy the BA equations 
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Where i = 1, . . . , N, a = 1, . . . , Mi and the parameters 
|Aa^| with I = 1,2, ...,k — 1 are the spin rapidities 

while we denote Xa^ = k a , A^' = and d = c/2. The 
BA quantum numbers Mj = U ~ * + l)iVj+i with 

M K — in the above BA equations ©. 



A. Repulsive regime 

The fundamental physics of the model are determined 
by the set of transcendental BA equations ^ which can 
be transformed to generalised Fredholm equations in the 
thermodynamic limit. The Fredholm equations for re- 
pulsive and attractive regimes are significantly different. 
From the BA equations the quasimomenta {ki} are 

real, but all j A^ j are real only for the ground state. For 
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where 1 < m < n — 1 and /3q = l/(27r), ro(fe) is the 
particle quasimomentum distribution function whereas 
r m (k) with m > 1 are the distribution functions for the 
k — 1 spin rapidities. The kernel function Kg(k) is defined 
as 



K e (k) = 



1 



Ic 



2tt {£c/2y 



(5) 



here c > for repulsive interaction and c < for at- 
tractive interaction. Follo wing the method used for the 
two-component Fermi gas [lq . we rewrite the Fredholm 
equations Q as 



m— 1 

(fc) = A)-5Z / *T m - s (fc-A)r s (A)dA 

s=0 ./|A|>B. 



A-!(fc-A)r m+1 (A)dA, 



(6) 



where < m < K— 1. The associating integration bound- 
aries -B m are determined by the conditions 



m e = = / r e (k)dk, 0<£<k-1, (7) 

L J-B e 

where Mo = N is the total number of fermions, N l = 
Me— i — Mg is the number of fermions in the £-th hyperfine 
state. Here M K = 0. The ground state energy E per unit 
length is given by 



E= [ k 2 r (k)dk. 

J-Bo 



(8) 



The model has SU(k) symmetry in spin sector and U(l) 
symmetry in charge sector. Therefore, the quantum num- 
bers of each spin states are conserved. Thus the sys- 
tem has k chemical potentials {fie} in regard to these 
conserved numbers N l . The ground state energy (JSJ is 
a smooth function of the densities of n = N e /L with 
I = 1,2, ... ,k for unbalanced case. In the grand canoni- 
cal ensemble, we can also get the chemical potentials fie 
via m — dE/dn l . 
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B. Attractive regime 

In the attractive regime, it is found that complex string 
solutions of kj also satisfy the B A equations 0, @, [13] ■ 
The quasimomenta kj may appear as bound states of 
m-atom up to length 2, . . . , k. A bound state in quasi- 
momentum space of length m takes on the form 

jgj = A (m-i) + l[m + l _ 2j) | c /| + Qicxpi-SL)), (9) 



where j = 1, . . . , m. The number of bound states with 
length 1 < m < n is denoted as N m . Its real part is 
A^™ 1 1 \ A k a bound state of m-atom will be accompanied 
by a Aq ■* string of length m — 1, a Xa^ string of length 
771 — 2 and so on until a Ai™ ^ string of length 1. Each 
accompanying string state in A^-space, A^-space, . . . , 
/ \(m- 1 )_ S p ace w in share the same real part AcT* , see 
[S H G3 and the second reference in [l2| . The unpaired 
atoms have real quasimomenta fcj's. 

From these quasimomentum bound states, the Fred- 
holm equations for the model with an attractive interac- 
tion are given by d, @ 



p m (X) = mf3 + £ £ / Ks 

(A-A) Ps (A)dA 

r=l s=r J — Qs 

+ K 8 . m (X-A)p s (A)dA, (10) 

s=m+l J~Qs 

where p± (k) is the density distribution function of single 
fermions, whereas p m (k) is the density distribution func- 
tion for the bound state of m-atom with 1 < m < k. The 
total number of fermions is given by N = X)m=i m ^m- 
The integration boundaries Q m , characterizing the Fermi 
points in each Fermi sea, are determined by 



—p- = / p m {k)dk. 
L J-Q m 



(11) 



The ground state energy per unit length is given by 



m(m 2 — 1) 
12" 



c 2 ) Pm (k)dk. (12) 



In the attractive regime, for convenience, we can de- 
fine the effective chemical potentials for the cluster bound 
sates p m = p+H m /m+(m 2 — l)c 2 /12 with m = 1, . . . ,k, 
where H m is the effective magnetic field (or Zeeman split- 
ting parameter) for the bound state of m-atom, see [12j . 
The Zeeman energy per unit length can be written as 
E z = 2r?i=i HmN m /L. In this regime, the effective 
chemical potentials for the bound states of different sizes 
can be derived from the energies of the ground state, i.e. 



1 d 

to dn„ 
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sk 2 p s {k)dk. 



(13) 



Following [T2| , the full phase diagrams of the model can 
be determined by the field-energy transfer relations 

H m = — m{n 2 - m 2 )c 2 + m (p m - Ha) (14) 

with m = 1, 2, . . . , k. In the next section, we shall de- 
rive the explicit form of the ground state energy that 
can give highly accurate effective chemical potentials p m 
to determine magnetism and full phase diagrams of the 
model. 

Similarly, the Fredholm equations can be rewritten as 
p m (\) = Po- [ K 1 (\-A) Pm - 1 (A)dA 

J\A\>Q m -i 

K rQ s 

+ Yl / K s - m {X-A) Ps (A)dA. (15) 

We see that there is a particular mapping between the 
two sets of Fredholm equations © and (fT5|) . i.e. 



— > 



C ->■ 



(16) 



This connection is useful in the analysis of the symmetric 
structure between the two sides. However, we find that 
the Fredholm equations ([5]) for the repulsive regime and 
(|15[) for the attractive regime do not preserve the map- 
ping which exists between the two sets of the Fredholm 
equations for the two-component Fermi gas 16]. This 
leads to particular intricacies in the analytical behaviour 
of the ground state energy as the interaction strength 
vanishes. In the present paper, we mainly concentrate 
on the solutions of the two sets of the Fredholm equa- 
tions for multicomponent Fermi gas with attractive and 
repulsive interactions. 



III. ASYMPTOTIC SOLUTIONS OF THE 
FREDHOLM EQUATIONS 

A. Strong repulsion 

For the balanced case, i.e. M m — M m+ i = N/k, the 
following theorem shows that all integration boundaries 
B m —s- oo except for m = 0, namely, there are no finite 
Fermi points (without chemical biases between different 
spin states). 

Theorem: For the K-component fermion system with 
the repulsive delta-function interaction, B m — > oo for 
m = l, — 1 if the relation 



M TO _i - M„ 

holds, and vice versa. 



Mr, 



M, 



m+l 



(17) 
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Proof: Integrating both sides of the Fredholm equations 
with infinite boundary (where m ^ 0), we obtain 



r m {k)dk 



M m M m +i 
L L 



If the condition (JTTJ) holds, thus we have 



r m (k)dk 



L 



(18) 



(19) 



where the integration boundary B rn = oo is inferred. 
Conversely, if B m with m > 1 tend to infinite, from Q 
and ([7]) we have 



Mm 

L 



r m (k)dk 



M m _i 



L 



M, 



m+l 
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For the balanced case, the integration boundaries B m 
with m > 1 are infinitely large. The fermi momentum 
-Bo is always finite. For convenience, we introduce the 
notation ro(fc) = ro in (k) + fo out (fc) where 



roi„(fc) = 

^Oout(fc) = 



ro(A:) 




ro(*) 



for |fc| < B 

for fcj > B 

for |fc| < B 

for \k\ > B 



Taking the Fourier transformation with the Fredholm 
equations (Q} for the balanced case, we obtain the fol- 
lowing relations 



-Fm(w) 



£ 

s=0 



-(re-m+2s)c|w|/2 



(20) 



where 2 < m < k — 1 and F re = 1. We denote the 
Fourier transform F[r m {k)\ = f m {uj). From the Fred- 
holm equation ([3]), we see that a closed form of the dis- 
tribution function r\{k) is essential for the calculation of 
the ground state energy. After a straightforward calcu- 
lation with the relations (|2"0|) . we obtain the closed form 



fo in (w) sinh [|(k — m)|w|c] 



sinh [i«|o;|c] 
that gives the distribution function 



Zir 



? m (uj) e - iulX duj. 



(21) 



(22) 



In the above equations f 0in (u>) can be determined from 
the following relation 



7(A) 



Ki{k - X)r Q (X)d\ 



-Bn 



Ki(k- \)r 0in (X)d\, 



(23) 
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FIG. 1: The ground state energy per length vs 7 = cL/N in 
natural units of 2m = ft = 1 for the balanced multicomponent 
Fermi gas: n — 2, 4, 10 and for spinless Bose gas in the whole 
interacting regime. For the repulsive regime, the discrepancy 
between the energies of the k — 10 Fermi gas and the spinless 
Bose gas is negligible. As k -> 00, the ground state energy of 
the balanced repulsive Fermi gas fully coincides with that of 
the spinless Bose gas [l9l ]. 



where their Fourier transform reads 7(w) = 
e~ c !"!/ 2 fo in (lj). The ground state energy per unit 
length ([5]) can be calculated either numerically or 
asymptotically. In Fig. [TJ we plot the ground state 
energy of ID balanced multicomponent Fermi gas for 
k = 2, 4, 10 and for spinless Bose gas [2(|. We see that 
as k —¥ 00, the ground state energy of the balanced 
K-component gas with repulsion coincides with that of 
the ID spinless Bose gas [l^. In order to capture the 
statistical nature for this connection, we shall calculate 
the first few terms of the expansion of the ground state 
energy. 

The strong repulsion condition, i.e., cL/N 3> 1, natu- 
rally gives c 3> -Bo 5 where -Bo is proportional to the Fermi 
velocity kp = nir. Following the method used in (l6| . we 
take a Taylor expansion with respect to A for the kernel 
function K\{k — A) in cq. (|2"3"|) . Using the relations ([7]), 
©, (22J and (E2D, we find 



fi(u) 



Euj 2 



F 2 (uj) 
Fi(u) 



Euo 2 



(24) 



In the above equations, we used the following formulas 



c 2 /4 + fc 2 ) 2 ] = 2nc- 3 e- 2 ^ 1 [c\u\ + 2} , 

Tic < - <- ^ - + %c\u)\ + 12] 



f (c 2 /A + k 2 y 3 -■ ' — 1 1 



Substituting (|24|) into the Fredholm equation ([3]), we 
obtain an asymptotic form of the density distribution 
function 



ro(*) 



1 

2^ 



jvy (fc) Br 2 (fc) 



2tt7 



4tt 



+ 0(c" 4 ), (25) 
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where the function 



Y a (k) 



u a F 2 (u) 



du. 



(26) 



After some algebra, we obtain the two functions used in 
2Z X 2Z 3 k 2 , nt _ 4 . 



r (fc) = — -— ^ + o(c" 4 ) 



c 
4Z 3 



with 



Z 3 



1 

K 

.-3 



+ 0(c- 4 ), 



^(-) + C 



C(3,-)-C(3) 

K 



(27) 



Here C( z )<?) an d CM are the Riemann zeta functions, 
tp(p) denotes the Euler psi function, C denotes the Eu- 
ler constant. When k = 2 we have Z\ = In 2, Z% = 
(3/4) C (3) that are consistent with the result given in (16| . 

Substituting (|2"5)l into © and ([5]) and solving by it- 
eration, the Fermi boundary Bq and the ground state 
energy of balanced ^-component Fermi gas with a strong 
repulsion are given explicitly (up to order 0(c~ 3 )) 
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2Z X 
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AZx 
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4Z 3 7T 2 
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Z\- 



Z 3 7T 2 

15 



(28) 



where the dimensionless interaction strength 7 = c/n. 
The ground state energy (f2"5)) with k = 2 reduces to the 
result given in [16]. In view of the quasi-momentum dis- 
tribution ro(fe) (US]), the system can be taken as an ideal 
gas with a less exclusive fractional statistics than the 
Fermi statistics. It is interesting to see that the ground 
state energy (j2"5|) reduces to the energy of the spinless 
Bose gas (see (10) in [l8[) as n —> 00, also see [lj|. Here 
we find that lim K _ i . 00 Z\ = IuHk-^oo Z3 = 1. A signifi- 
cant interpolation between Fermions and Bosons can be 
conceived from the parameters Z\ and Z 2 in (|2^|) that en- 
codes the quantum statistical and dynamical effects. In 
Fig. [2] we see a good agreement between the asymptotic 
expansion result (|29p and numerical result obtained from 
the Fredholm equation ([3]) with the density distribution 

In general, the boundaries B m with m > 1 are hard 
to be estimated for the imbalanced case in a strong re- 
pulsive regime. However, for strong repulsion, the spin 
sector is dominated by the density-density interaction. 
Therefore, the spin polarization is not essential. We con- 
sider the high polarization case where M m _i — M m <C N. 
In this case, the condition c 3> B m always holds. Thus it 
is straightforward to calculate density distributions from 
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FIG. 2: The ground state energy per length vs 7 = cL/N in 
natural units of 2m = h = 1 for the balanced multicomponent 
Fermi gas: k = 2, 4, 10 and for spinless Bose gas. The solid 
lines are the numerical solution obtained from the two sets 
of Fredholm equations. The dashed lines are plotted directly 
from the asymptotic expansion ground state energies (I29|) . 
(|44[) for a repulsive interaction and (|55[) . (159|l for an attrac- 
tive interaction. High precision of the asymptotic expansion 
ground state energy is seen through out strongly and weakly 
coupled regimes. 



Q and Q with a proper Taylor expansion (up to order 
0(c- 3 )),i.e. 



ro(fc) 
ri(fc) 

r m (k) 



4P 
„2 



1 2Mi 

27T C7Ti 

— - \2N - Mi + 2M 2 ] 

[87V - Mi + 8M 2 ] - 



—T [2M m _i 

[8M m 



c 3 irL 



M m + 2M m+ i] 
1 - M m + 8M m+ i] 



(30) 



where m = 2, . . . , k. It is also straightforward to cal- 
culate the Fermi boundary B and E from ([7]) and (jHJ), 
i.e. 



Bo 



n7r < 1 



4mi 

c 



16to 2 



64mf 



167T mirr 



3c 3 



1 



+ 0(c" 4 ), 



8mi 48m 2 256mf 



3 L 

327T 2 min 2 
5c^ 



+ 0(c^ 4 ), 



(31) 



where mi = M x /i with Mi = J2j=l Ni+1 - For thc 
highly polarized case, Mi <C TV, the ground state energy 
E per unit length, up to 0(c~ 3 ), solely depends on the 
quantum number Mi . The spin states are not essential in 
this strongly repulsive regime due to the freezing of spin 
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transportation. This result agrees with the energy of the 
two-component Fermi gas with a strong repulsion, see 
[lH . It clearly indicates that strong repulsion suppresses 
the spin effect. 



B. Weak repulsion 

For the balanced case with weak repulsion, i.e., 
cL/N <C 1, all B m with m > 1 tend to infinity, so the 
Fredholm equations © become 



r (fc) =$o+ K x {k - X)n{X)dX, 

J — oo 

/oo 
K m (k - X)r 0oM (X)d\ 
-oo 

POO 

+ Kx{k - X)r m+1 (X)d\, (32) 



where 1 < m < k — 1. By iteration with r\, r 2l ■ ■ ■ , 
the first equation of (I3"2l becomes 



ro 



K 1 /.OO 

(k) = Kpo~Y, K 2s (k - X)r 0o nt(X)dX. (33) 

„_ i J —oo 



Using a Fourier transformation, we easily prove 



r o „t( fc ) = A) 
i?(fc) 



R(k - A)r 0in (A)dA, 



-So 



1 

2^ 



dw. (34) 



Here we see that R(k) = 0(c). Substituting the leading 
term of rQ axxt (k) — /3q into (|33[) , we further obtain 



ro(*) 



K-l r 

--— W 

2tt 2tt 2 ^ 1 

s=l k 



arctan 



arctan 



B + k 



B -k / 
+ 0(c 2 ) (35) 



for the region \k\ < Bo. The energy can be calculated 
from the equations (|7|) and © with the distribution func- 
tion (E3 



E 



ir 2 n 3 



3k 2 



+ c(k- l)n 2 /n + 0(c 2 ). 



(36) 



When k = 2 the energy (|3"6")) reduces to the result given 
in This result clearly indicates a mean field ef- 

fect among the balanced K-component weakly interact- 
ing fermions. The kinetic energy part in Q36p vanishes as 
k — y oo. The energy (|36[) as k — > oo is the same as for 
spinlcss Bosons with a weak repulsion. 

For the imbalanced case with weak repulsion, we as- 
sume B m _i > B m with 1 < m < k — 1. The calculation 
of the ground state energy is very complicated because 



the k integral equations are coupled with each other. We 
have to separate the integration intervals case by case. 
Using the Fredholm equations ([5]), we calculate the fol- 
lowing integral in different regions. We denote the in- 
tegral I m = J^g r m (X)Ki(k — X)dX. For the region 
|fc| > B m , we find 



£cB n 



27T 2 (fc 2 - B 2 n 



E 

s— m+ 



cB x 



i 27T 2 (fc 2 -B 2 S ) 



After a lengthy calculation, we calculate the integral for 
the region |fc| < B m 



In 



£cB n 



2^ 2 {B 2 m -k 2 ) 



E 



(£-m + s + l)cB s 

2tT 2 (B 2 -k 2 ) 



r m+1 (X)K e+1 (k - X)dX + 0(c 2 ). (38) 



For Bp < \k\ < Bp-i with p > m, we further calculate 



(p — m)(£ + ,s — m + l)oE> s 



2tt 2 (B 2 



E 



[{p -s)(£ + s-m + l)- l]cB s 



(£ + p — m)cB p 



2n 2 {B 2 -k 2 ) 

K-l 

E 



cB s 



2^ 2 (fc 2 -B 2 ) s +^ 2 (fc 2 - B 2 ) 



(39) 
+ 0(c 2 ). 



With the help of these formulas, we are able to evaluate 
the order of r m {k) in the Fredholm equations (JBJ) 



r m {k) = 0(c) 
r m (k) = fa + 0{c) 



for fc| > B m -\ 
forB TO+ i < \k\ < B m -i 



r m (k) = (p- m)j3 + 0(c) for 



B p < \k\ < B p _i 



(.40) 



p > m + 1 

From ([7]), we obtain the integration boundaries via 



Mm 

L 



M, 



7T L 

+ E ln 



m+l 



B s — B„ 



B m — B s 



B m + Be 



B s + B„ 
+ 0(c 2 ), 



(41) 



s— m+l 

where < m < k — 1. In order to calculate the ground 
state energy, we also need a lengthy calculation of the 
integral 

"m— 1 



X z r m (X)dX = 



Bt cB n 



E ». 



s— m+l 



3tt tt 2 

( m— 1 



-±-\Y B 2 S In 



E B > 



s— m+l 
r-B m+ l 



27T , 

k s=0 

B m — B 



. s=0 

B s — B n 



B s + B m 



B m + B s 



X 2 r m +i(X)dX + 0(c 2 ). 



(42) 
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Using ([8j and (|42)) . the ground state energy i£ per unit 
length is given by 



E 



E 

m=0 



B 3 

rn 

3tt 



K — 2 K. — 1 



(Bt 



^E E [4^3 

m— r— m+1 

£>„,. — -B 



B„ In 



+ 0(c 2 ). (43) 



Substituting (j4"Tj) into (|43|) . we finally obtain the ground 
state energy of the K-componcnt Fermi gas with weak 
repulsion 



^ K — 1 

E = gX^- 1 

1=1 



3^2 



rrii) 7r 



(44) 



+ 2C E E ( m i-l - m i)( m j-l - m ]) + 0{C 2 ), 
i=l j=i+l 



with mo = n and m K = 0. Here the linear density 
ri = N/L and the quantum numbers m, = Mi/L. If we 
introduce polarization = N l / L with i = 1,2, 1, 
where iV* is the number of fermions in the ith level. Thus 
we have a simple form of the ground state energy per 
length 



i=i 



2c E E P*P 3 +0(c 2 ). (45) 



i=l 



The first part is the kinetic energy of the K-component 
fermions whereas the second parts is the interaction en- 
ergy. This result is valid for arbitrary spin imbalance 
in the weakly repulsive regime. This result is in a good 
agreement with the numerical calculation, see Fig. [2] 
For the balanced case, i.e. M m — M rn+ i = N/n, (|4"4")l 
reduces to the energy (j3"6| . It is interesting to note that 
the ground state energy (|45"1) presents a mean field theory 
of the two-body s-wave scattering physics. 



C. Weak attraction 

In weakly attractive coupling regime, the Fredholm 
equations give the distribution functions of clusters of 
different sizes. For weak attraction, i.e., \c\L/N <C 1, the 
two sets of the Fredholm equations for repulsive and at- 
tractive regimes preserve the symmetry (|16p . Therefore, 
the calculation of the ground state energy E per unit 
length for weak attraction is similar to that for weak re- 
pulsion. For the balanced case with a weak attraction, 
N m = for m = 1, ...,« — 1 and N K = N/k. Thus, we 
have the condition Q m = for all Fermi points except 
Q K . In this case, the ground state is a spin singlet state. 
The Fredholm equations (|10p for the spin neutral bound 
state of a K-atom becomes 

Pk (\) = k(3 q + J2 K 2s (X-A)p K (A)dA. (46) 



s=l 



Under the mapping (|16[) . the Fredhlom equations (|46[) 
for weak attraction map to (|33p for weak repulsion. In 
order to estimate the contribution of density distribution 
p Kia , we take the Fourier transformation of (|46[) and then 
prove that 



T(A) 









/ T(X 










2^ J 


— oo 



E« 



-s|ccj| 



duj, 



where T(X) = 0{c). For |A| < Q K , we thus have p K (X) — 
p Kin (A) = /9o + 0(c). Substituting this leading order into 
(|46l) . we have 



p K .(X)dX = ^ + 
s\c\ 



(k-1)\c\ 



K-l 

V 

^ 2ktt 2 

S— 1 



In 



7T 2tt 2 
s 2 \c\ 2 +4Ql 



+ 0(c 2 ). (47) 



Similarly, from (IT2"T) , we calculate the ground state energy 
of the balanced gas 



E = 



|c|( K -l)n 2 



3k 2 



0{c 2 



(48) 



We see that the energy of the balanced gas with weak at- 
traction (|48p and with weak repulsion (|36p continuously 
connect at c — > 0. It is also seen that (|4"5)l with k = 2 is 
consistent with the result given in (l6| . 

For the imbalanced case, the ground state has clus- 
ter states of different sizes. Therefore we assume an 
ansatz in order to calculate the integration boundaries 
for each of the Fermi seas, i.e. Q m > Q m +i with 
1 < m < k — 1. Similarly, we should calculate the inte- 
grals I m = £ J Q Q m ( Z^2)i n +(l d -kr for different regions. 
After some algebra, we find 



En. 



£\c\Q r 



\c\Qs 



2ir 2 (X 2 -Q m ) s= ^ +i 27r 2 (A 2 -Q 2 ) 
for the region |A| > Q m and 

T _ „ Z\c\Q m t\c\Q m -\ 

J-rn — PO 



0(c 2 )(49) 



2^(Ql -A 2 ) 2n 2 (Q 2 m _ 1 



A 2 ) 



+ 



ro— 2 



\c\Q s 



^ 2vr 2 (Q 2 -A 2 ) 

■Qs 



(50) 



E 

s—m+1 



K e+s - m (X - A)p s (A)dA + 0{c 2 ). 



for the region |A| < Q m . In the region Q m +i < |A| < Q m , 
we further obtain 



m-2 



Im = A) + E 



\c\Q s 



^ 2^ 2 (Q 2 -A 2 ) 2 7 r 2 (Q 2 l _ 1 -A 2 ) 
t\c\Q m ^ (£ + l)\c\Q s , ^ 2 



2 7 r 2 (Q 2 l -A 2 ) s £^ +1 2vr 2 (A 2 - Q 2 ) 



0(ci51) 
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For |A| < Qm+i, we find 

(e + i)\c\Q m+1 



{2i+l)\c\Q T 



2v 2 (Q 2 m+1 -\ 2 ) 2^{Q 2 m -\ 2 ) 



{l-l)\c\Q m -i 



2ir 2 (Q 



m— 1 



A 2 ) 



+ 0(c 2 ). 



(52) 



From these equations (|49ll52p . we are able to evaluate 
the leading order contributions in the density distribution 
functions, i.e., p m (X) = (m—p+ l)/?o + 0(c) for Qp+i < 
|A| < Q p with p < m and p m (A) = O(c) for |A| > Q m . 
From equation we obtain the integration boundaries 



Nrn Qrn^ A _ Jcj_ j ^ 

s— m+1 k s — 1 



Qr 



E 

s— m+1 



+ 0(c 2 ). 



(53) 



From (fT2")l . the ground state energy E per unit length is 
given by 



E 



E 

m — 1 



3vr 



2|c| 



^ ] ^ ' QsQr 



(54) 



s— 1 r— s+1 



\c\ 
2tt 2 



53 E Wr + Q 2 )ln 



S — l T — S + l 



0(c 2 



Substituting the integration boundaries Q m from (|53[) 
into (|54l) . we thus obtain the ground state energy of the 
K-component Fermi gas with a weak attraction 



£ = ^E^ 2 -2|ciE E 



PiPj, 



(55) 



i=l j=i+l 



where we introduced polarization p.; = N 1 /L with i = 
1,2, ...,«;. The first part is the kinetic energy of the 
K-component fermions whereas the second part is the 
interaction energy. This result is in a good agreement 
with numerical calculation, see FigJ5] For the balanced 
case, i.e. TV 1 = N/k, (|55[) reduces to the energy (|4"5| . 
From the energies (|45]) and (|55p . we see that the ground 
state energy of the K-component gas with arbitrary polar- 
ization continuously connects at c = 0, precisely speak- 
ing, the ground state energy per unit length and its first 
derivative is continuously connected at c = 0. Again, the 
ground state energy (|55[) presents a mean field theory of 
the two-body s-wave scattering physics in weak interact- 
ing regimes. 



D. Strong Attraction 

Universal low temperature behaviours of isospin S — 
1/2, 1, 3/2, . . ., (k— l)/2 interacting fermions with an at- 
tractive interaction in ID shed light on the nature of 



trions, pairing and quantum phase transitions. The ex- 
istence of these internal degrees of freedom gives rise 
to some exotic superfluid phases. These models exhibit 
new quantum phases of matter which are characterized 
by bound states of different sizes underlying the symme- 
tries. For example, three-component ultracold fermions 
give rise to a phase transition from a state of trions into 
the BCS pairing state under external fields. Recently, 
considerable interest has been paid to the low dimen- 
sional strongly interacting fermionic atoms with high spin 
symmetries. 

For a strong attraction, i.e. \c\L/N 1, the bound 
states of different sizes form tightly bound molecules of 
different sizes with binding energies £^ = £(£ 2 — l)c 2 /12, 
where I = 2, 3, . . . , k. In this regime, all the Fermi mo- 
menta of the molecules are finite, i.e. |c| 3> Q m with 
m = 1, . . . , K. Here Q\ charaterizes the Fermi momentum 
of the single spin-alined atoms. In the canonical ensem- 
ble, external fields or nonlinear Zeeman splittings trigger 
rich quantum phases and magnetism [l2| . A closed form 
of the ground state energy with polarization is essen- 
tial to work out phase diagrams and magnetism at zero 
temperature. In view of the strong attraction condition 
|c| 3> Q m and following the method [l6j], proper Taylor 
expansions can be carried out for the three supplemen- 
tary formulas 

^(A-A).A- 4Q 16 W + 3 * 2 ) 



n£\c\ 

64Q(Q 4 + 10Q 2 A 2 + 5A 4 ) 

5n£ 5 \c\ 5 
2 

(A - A)K e {\ - A)dA 

Q 

4QA 16QA(Q 2 + A 2 ) 
~wi\c\ + ir£ 3 \c\ 3 

A 2 K e (X — A)dA 

Q 

4Q 3 16Q 3 (3Q 2 + 5A 2 ) 



3tt£ 3 | 



+ 0(c" 7 ), 



+ 0(c-% 



+ 0(c~ 5 ). 



In the above calculation, we denote the kernel function 

*w - 1 f|c| 



3n£\c\ 15tt£ 3 |c| 3 
we dt 

, . . The following calculations are 
2tt (£c/2) 2 + k 2 b 

valid for arbitrary polarization. For a strong attraction, 
all terms in the Taylor expansion of the Fredholm equa- 
tions are converged well, see the method proposed in [16| . 
From the Fredholm equations (|10[) and the boundary con- 
ditions (TlTl) , it is straightforward to obtain the following 
expression (up to the order 0(c~ 3 )) 

16tt 2 



Nm 

L 



mQ r , 



1 - 



mL\c\ 



:F„ 



Q„ 



A 2 p m (A)dA 
Wtt 2 



3tt 



3to 3 j L 3 |c| 3Gi,; ( 
4 



mlc 



rF„ 



15m 3 L 3 |c| 3 



(56) 
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where we denoted the functions 



771 — 1 S 

F m = YY- 

s — l r—1 r—1 



m-1 N 
Or 



k m 



E 



^— ' 2r + s — to — 2 ' 

s— m+1 r—1 



2. 2. s 2(2r-s + TO-2) 3 ' ^ (2r) 3 
s 2 ^iv s + TO 2 7V 3 



k m 



E 



v 



' , ^ s 2 (2r + s - to - 2) 3 ' 

+1 r—1 v / 



^^9s 2 N^N s + 5m 2 N 3 



2.2. 



14N 3 



^ ^ s 2 (2r-s + TO-2) 3 ^ (2r) 3 
9s 2 iV 2 iV s + 5TO 2 7V 3 



k m 

E 

s— m 



•> E s-(2r * - ■ m 2i :i 

+lr=l v ' 



Here Ni with i = 1,2..., are the numbers of the cluster- 
state of j-atom. 

After a lengthy calculation, we obtain explicit forms of 
the Fermi momenta and the energies of the molecules of 
different sizes for a strong attraction 



N„ 



1 



mL [ toL|c 
16 



- F 



16 



-F 



3m 3 L 3 |c| 3 
n 2 N 3 



1 



3toL 3 1 mLlc 



m 2 L 2 \c\ 2 "" 
[UFl - G m n 2 ] | , (57) 
s 

-F 



256 



L 2 |c| 2 m to 3 L 3 |c 



3 m 



16tt 2 



i 3 L 3 | 



[-G m + G m /15]|. 



(58) 



Thus the ground state energy of the gas with arbitrary 
polarization per length in strong attractive regime is 
given by 



(59) 



where the binding energy of the molecule state of £-atom 
is given by e[ e) = l{l 2 - l)c 2 /12. For k = 2, it cov- 
ers the result obtained for the two-component Fermi gas 
given in (l6| . From the discrete BA equations, one of 
the author and coworkers [12j derived the ground state 
energy of K-component strongly attractive Fermi gas for 
up to the order 0(l/c 2 ). Here we obtained a more ac- 
curate ground state energy ([55)1 of the K-component gas 
with arbitrary polarization from the analytical study of 
the Fredholm equations. We noticed that the ground 
state of strongly attractive K-component Fermi gases 
can be effectively described by a super Tonks-Girardeau 
gaslike state via a proper mapping 21]. The explicit 



form of the ground state energy can be used to study 
magnetism and full phase diagrams of the model from 
the relations (ITBl and in a straight forward way. 
The effective chemical potentials may be calculated from 
Mm = d Ei) I (mdn m ), where n m = N rn /L. From 

these effective fields, the phase diagrams can be deter- 
mined by the energy- field transfer relations (I14[) . see re- 
cent study of the three-component attractive Fermi gas 

MM . 



IV. LOCAL PAIR CORRELATIONS 

There has been a considerable interest in studying uni- 
versal nature of interacting fermions. Remarkably, Tan 
[23j showed that the momentum distribution exhibits 
universal C/k decay as the momentum tends to infin- 
ity. Here the constant C is called universal contact that 
measures the probability of two fermions with opposite 
spin at the same position. For ID two-component Fermi 
gas 24[, the universal contact is obtained by calculat- 
ing the change of the interacting energy with respect to 
interaction strength by Hellman-Feynman theorem, i.e. 
C = ^^n-|-n^<?j 2 j(0). The local pair correlation are 
accessible via exact Bethe ansatz solution through the 
relation g| 2 |(0) = dE/dc. Here E is the ground 

state energy per length. 

In a similar way, for ID K-component Fermi gas, there 
exists a ID analog of the Tan adiabatic theorem where 
the universal contact is given by the local pair correla- 
tions for two fermions with different spin states. The 
two-body local pair correlation function is similar to the 
calculation of the expectation value of the four-operator 
term in the second quantized Hamiltonian. For a ho- 
mogenous and balanced K-component Fermi gas, the lo- 
cal pair correlation function is given by 



dE 



(k — l)n 2 dc 



(60) 



with k > 1. Where E is the ground state energy per 
length. From the asymptotic expansion result of the 
ground state energy of the balanced Fermi gas obtained in 
the above section, we easily find g y a ' a ,(0) — »• 1 as \c\ — > 0. 
From the ground state energy (|2"9")l of the balanced gas 
with a strong repulsion, we have the local pair correlation 



4k7T 



3(k- 1)7 2 



Zi 



6Z 2 

1 



24 
— 2 



Zf 



15 



2 \ 1 



(61) 



This local pair correlation reduces to the one for the spin- 
less Bose gas [H,[25| as k —> 00. 

For strong attractive interaction, we obtain the ground 
state energy of the balanced gas from the result 



E = 



3k 4 
32B 



1 



4A K 
* 2 l7l 



12A 2 



K 4 7 2 



32^ 

K 6 |7| 3 



15k 6 7 3 



- n K .e 



(«) 



■0(l/ 7 4 ), 



(62) 
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FIG. 3: The local pair correlation <?' 2 , (0) vs 7 = cL/N 
in natural units for the balanced multicomponent Fermi gas: 
k = 2,4, 10 and for spinless Bose gas. The solid lines are 
the numerical solutions obtained from the two sets of Fred- 
holm equations. The dashed lines are the asymptotic limits 
obtained from the first term in (|63|) . All local correlations 
pass the point g^,(0) = 1 at vanishing interaction strength. 
However, the local pair correlation diverges for the attractive 
Bose gas. 



where A K = J2r=l V r an d B K = J2r=l V 7 " 3 - From the 
relation we obtained the local pair correlation for 

two fermions with different spin states in a strong attrac- 
tive regime 



7( 2 ) 



«(« + i)M 

6 
24 



4^ 



3ft 5 (ft 
B, 



6AI 



A 3 - — 
k 4 7 2 V " 15 



1)7 2 \ k2 \i\ 
+ 0(l/ 7 4 ). (63) 



We see that the local pair correlation becomes divergent 
as k — )■ 00, i.e. goes to the limit for the attractive bosons. 
In Fig(3j we present the numerical solution of the lo- 
cal pair correlation for two fermions with different spin 
states in the multicomponent Fermi gas. In the whole 
interacting regime, the local pair correlations for the bal- 
anced K-component Fermi gas tends to the limit value for 
the spinless bosons as k — > 00. High precision of these 
asymptotic expansion local pair correlations is seen in 
strong and week interaction regions, see Figf4] 



V. CONCLUSION 

We have analytically studied the Fredholm equations 
for the ID fc-component fermions with repulsive and at- 
tractive ^-function interactions in four regimes: a) strong 
repulsion; B) weak repulsion; C) weak attraction and D) 
strong attraction. Solving two sets of the Fredholm equa- 
tions (01 and (|10p . we have obtained the first few terms 
of the asymptotic expansion for the density distribution 
functions and the ground state energy E per unit length 



FIG. 4: The local pair correlation g ,(0) vs 7 = cL/N 
in natural units for the balanced multicomponent Fermi gas: 
ft = 2, 4, 10 and for spinless Bose gas. The solid lines are the 
numerical solutions obtained from the two sets of Fredholm 
equations. The dashed lines are the asymptotic expansion 
result obtained from (1611) and (16311, 



for both balanced and unbalanced cases in these regimes. 
We summarize our main result as following. 

A) For the strong repulsive regime, the ground state 
energy of the balanced ft-component gas has been given 
in ([29]) up to the order of 1/c 3 . It presents a universal 
structure in terms of the parameters Z\ and Z3 (|27|) that 
characterize quantum statistical and dynamical effects of 
the model. It also confirms Yang and You's result [l9| 
that the energy per particle as k — > 00 is the same as for 
spinless Bosons. In this regime, we have also obtained 
the ground state energy for the highly imbalanced case 
(|3"1"|) . It clearly indicates that the internal spin effect is 
strongly suppressed in this regime, i.e. the energy only 
depends on the quantum number M±. 

B) For the weak repulsive regime, we have calculated 
the ground state energy E (131)1) for the balanced case and 
(|44p for the imbalanced case. The ground state energy 
(l44l) presents a simple mean field theory of the two-body 
s-wave scattering physics. We see that the energy ([56]) 
for the balanced ft-component Fermi gas is the same as 
that for the spinless Bose gas as k — s> 00, consistent with 
the result in [19j . 

C) For the weak attractive regime, we have also calcu- 
lated the ground state energy E (|48|l for the balanced case 
and (|55|) for the imbalanced case. We have found that the 
two sets of the Fredholm equations for the model with re- 
pulsive and attractive interactions preserves a particular 
mapping (ITB1) . But they do not preserve the density map- 
ping found for the two-component Fermi gas [la ]. Our 
result shows that the ground state energy of the model 
with arbitrary spin state is continuous at c — > 0, as is the 
first derivative. 

D) For the strong attractive regime, the highly ac- 
curate ground state energy (up to 0(c~ 3 )) of the ft- 
component with arbitrary spin polarization has been 
given in ([55]) . The ground state of the largest cluster 
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state can be effectively described by the gas like state 
in the super Tonks-Girardeau gas of hard core bosons. 
From the result of ground state energy of the strongly 
attractive fcrmions, we can study magnetism and iden- 
tify full phase diagrams of the gas with SU(n) symmetry 
driven by the external fields. 

Furthermore, for the balanced case, we have obtained 
the first few terms of the asymptotic expansion for the 
local pair correlations of two fermion with different spin 
states in the four regimes. A numerical solution of the 
local pair correlation function has been obtained in the 
whole interacting regime. These results we obtained pro- 
vide insight into understanding quantum statistical and 
dynamical effects in ID interacting fermions with high 
spin symmetries. The explicit forms of the ground state 



energy of the multicomponent Fermi gas with polariza- 
tion give further applications to the study of quantum 
phase transitions, magnetism, and phase diagrams, which 
provides useful guides for future experiments with ultra- 
cold atomic Fermi gases with high symmetries. 
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